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FOREWORD 


This  report  presents  the  results  of  a  study  to  determine  the  magnitude  of  the  second- 
and  third-order  ionospheric  corrections  to  an  electromagnetic  wave  propagating  through  a 
spherically  symmetric  refracting  medium.  The  results  show  that  both  the  second-  and  third- 
order  correction  terms  must  be  included  in  a  study  of  higher  order  ionospheric  effects,  contrary 
to  previous  reports  that  the  second-order  term  is  negligible.  This  report  was  reviewed  by 
Dr.  Jeffrey  N.  Blanton,  Head,  Space  Flight  Sciences  Branch. 


Released  by: 

THOMAS  A.  CLARE,  Head 
Strategic  Systems  Department 
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INTRODUCTION 

The  phase  path,  L^,  is  the  curved  path  that  an  electromagnetic  wave  follows  when  traveling 
from  a  satellite  to  a  receiving  station.  The  path  is  curved  due  to  refraction  effects  of  the 
Earth’s  ionosphere. 

An  expansion  for  the  phase  path  length  takes  the  form 
=  p  +  L,  +  La  +  L3  +  ■  •  • 

where  p  is  the  vacuum  slant  range  from  the  station  to  satellite  and  L,,  La,  L3,  •••  first-, 
second-,  third-,  etc.  order  corrections  to  the  phase  path.  By  assuming  the  refracting  medium  to 
be  spherically  symmetric,  a  closed-form  integral  can  be  derived  for  L^  and  subsequently  for 
L,,  La.  L3,  etc.  The  clo.sed-form  integral  expressions  have  been  derived  by  Gibson'  for  the 
first-,  second-,  and  third-order  corrections  and  are  dependent  on  the  vacuum  zenith  angle  and 
the  index  of  refraction  model  used. 

The  first-order  correction  is  given  by' 


where  the  constant  k  is 

k  =  r  sin  / 

0 

with  r  the  distance  from  the  center  of  the  Earth  to  the  station  and  z  the  vacuum  zenith 
angle.  N(r)  is  the  index  of  refraction  of  the  ionosphere  or  plasma; 

'  1..  R.  (ahMM).  .SV)»/C  /  lor  an  l■.h'^■tromal;)u■lk■  Ui;rc  I'ropuKalma  Throufrh  a  Sphrricallv  Svninu  tric  Rffraclinn  .Mi\liiini.  Nav;il 

Surf.KC  WtMpoiis  (  etifor  fechnic.i>  Kepnrf  .NSVIX  /HI  I  K  .K)44  Dahlprcn.  VA  Junv.  1975. 
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n*(r) 


(2) 


where  cjp(r)  is  the  plasma  frequency  and^ 


(r)  = 


4ire^ 


m 


N(r) 


(3) 


N(r)  is  the  electron  density  in  electrons/cm^ ,  e  the  electronic  charge  in  statcoulombs,  m  the 
mass  of  the  electron  in  grams,  and  co  the  signal  frequency  in  radians/s. 

Substitution  of  Equation  2  into  Equation  1  shows  that  the  first-order  correction  goes  as 
1/cj^: 


L.  = 


roiHx)  rdr 

1/2 

J  V r^-k^ 


(4) 


SECOND-ORDER  CORRECTION 


The  second-order  correction  to  the  phase  path  length  accounts  for  the  effect  of  the 
Earth’s  magnetic  field  on  the  ionosphere. 

The  index  of  refraction  of  a  plasma  in  an  external  magnetic  field  is:^ 


n^(r) 


cjCco  + 


(5) 


where  cjg(r)  is  the  frequency  of  precession  of  a  charged  particle  in  a  magnetic  field: 


coalr) 


eB(r) 

me 


(6) 


D.  Jackson.  Classical  Electrodynamics ,  John  Wiley  and  Sons.  Inc..  New  York,  NY.  1975 
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In  Equation  6.  e  and  ni  are  as  defined  previously,  c  is  the  speed  of  light  in  cm/s.  and 
B(r)  is  the  component  of  the  Earth’s  magnetic  field  in  the  direction  of  signal  propagation. 
From  Equation  5  it  can  be  seen  that  electromagnetic  waves  of  different  polarizations  propagate 
differently. 


To  obtain  the  expression  for  the  second-order  correction,  expand  Equation  5; 


n^(r)  =  1  - 


1  ± 


co,,(r) 


CO 


n^(r)  -  1 


co„(r) 

1  .i. 

\ 

1  X 

CO 

- 1 

=  1  + 

cOjj/oj  +  co: 

-  coMr) 

P 

co^ (r)  CO 

p'^  '  H 

co^ 

"  C03 

4/CU3 


U) 


TIk'  kist  expression  substituted  into  E(|uation  I  modifies  the  first-order  correction; 


L, 


~C  -  co^(r)  rdr 

10  1  - "  " 

^  co^x/r'-k* 


■¥ 


f 


co^(r)  co„(r)  rdr 

P _ H _ 

cj3  \/r^-k^ 


1 


to*(r)  coj, (r) 
I’  H 


CO 


4 


rdr 


(7) 


file  first  term  in  tlie  e\|iansion  is  just  the  first-order  correction.  L, .  as  given  previously  in 
l  (|uation  4:  the  second  term,  which  gt)es  as  1/co^.  is  called  the  second-order  correction,  L2 : 
and  tlie  tliird  term,  which  goes  as  I/co'’,  is  a  part  of  the  third-order  correction  and  will  be 
discussed  later.  Note  that  the  secoiul-order  correction  is  tlie  first-order  correction  multiplied 
b\  tlie  precession  freciucncN  and  tliviiled  h\  the  oscillator  freiiuency: 


CO 


To  determine  tlie  magnitude  of  Li-  it  is  necessary  to  find  cOj^  ■ 


(8) 


I'o  find  tlie  precession  tre(|uene\.  an  exiuession  t’or  B(r)  is  needed  in  the  direction  of  signal 
propacalion.  p  (l  iciirc  1). 


3 
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FIGURE  1.  DIRECTION  OF  SIGNAL  PROPAGATION 

Given  the  Earth-fixed  coordinates  of  the  satellite  (x  ,  y^,  z  )  and  station  ( 
p  can  be  defined  by 

i?  -  (X  -  X  )  X  +  (y  -  V  )y  +  (z  -  z  )z 

S  (>  S  O  s  o 

Tlie  magnetic  field  of  the  earth  can  be  represented  in  geodetic  coordinates  as 
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is  the  component  of  B  in  the  direction  of  signal  propagation.  The  resulting  expression  for  is 

B  =  — —  Fb  l(x  -  X  )  cos  9  cos  0  +  (y  -  y  )  cos  0  sin  0  +  (/.  -  z  )  sin  0  ) 

IPI  L 

+  Br,  [-  (x  -  X  Isin  0  cos  0  -  (v  -  y  )  sin  0  sin  0  +  (z  -  z  )cos  0) 

+  B  [-  (X  -  X  )sin  0  +  (v  -  y  )cos  0| 

Also  reciuired  before  evaluation  of  Fquation  8  can  begin  is  a  model  for  the  electron 

density  profile,  N(r).  in  Fciuation  3. 


The  ionospliere  is  composed  of  several  distinct  regions  in  which  the  functional  form  of 
N(r)  varies  from  one  region  to  the  next.  Therefore,  it  is  necessary  to  separate  the  integral  in 
Hquation  8  into  a  series  of  integrals  over  the  different  ionospheric  regions: 


-TTC^  CJ. 


N|  (r)  rdr 


1:  Liver 


f  N|  (r)  rdr  r  N^. 

J  V'r-k^  J 

FI  Laver 


(r)  rdr 


f  — 

J 


(r)  rdr 


(r)  rdr 


Lower  F2  Layer 


mid-upper  1'2  Laser 


/N  (r)  rdr  /-N 

s/r^-k^  J  V  r 


N,,  (r)  rdr 


'  x/r'-k'  j  Vr'-k' 

1st  upjK’r  F2  2nd  upper  F'2 

exiionentially  exponentially 

decaying  layer  decaying  layer 


I  Vr^-k^ 

3rd  up|ier  F2 
exponentially 
decaying  layer 


For  the  numcric,il  analssis  jx’rformed  later  in  this  report,  the  functional  forms  for  tlic  corre¬ 
sponding  Nlrfs  were  taken  from  tlie  .Apidied  Research  Laboratories  (ARL)  ionospheric  model 
( Figure  2  ' 


1.  H.  (  hncJi  .iit.J  K.  \.  SJu-nhury.  ]>>rn<\rln  rw  AV'/l’f  f  rrnr  MmiiL  .Xpplii-iJ  Krsi-.jnl) 
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To  include  tlie  majinelic  field  interaction,  the  third  term  in  hquation  7  must  be  added  to 
the  above.  Substituting  for  n^(r)  -  1  and  adding  the  magnetic  field  interaction  term  leads  to 
the  following  expression  for  the  third-order  correction: 


I: 


rMr 

r  toMrirdrI 

iVr- 

J(r^-kM-’'*^ 

/J(r> 

rdr 


i 


-I 


cc'ir)  w:. (rl  rdr 

I’  u 


v''r'-k‘ 


(9) 


The  third-order  correction  goes  as  1  oj'* .  As  before,  each  integral  in  Ecjuation  9  must  be 
separated  into  a  series  of  integrals  over  tlic  different  layers  of  the  ionosphere. 


SINGLE-FREQUENCY  DOPPLER  DATA 


I’or  satellites  emitting  a  single  frequency,  fj. ,  the  first-,  second-,  and  third-order  iono¬ 
spheric  corrections  can  be  determined  directly  from  luiuations  4.  8.  and  9.  I'actoring  out  the 
transmitter  frei|ucncy  gives  tlic  following  convenient  form  for  the  phase  path; 

I  I  I 

Lp  -  E  +  L,  +  ^  L:  +  ”  Lj 
T  1  I 

for  Doppler  data,  the  time  deriv.itice  of  the  phase  path  is  needed  to  calculate  the  received 
frequenev.  f^ 


f  =  (1 

r 


elf, 


(10) 


where 


f- 


L,  + 


I  .. 
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Expunding  Equation  10  gives 


•  -  ■  ±  •  _L  •  J_  • 

-  c  '  cf^  ^  “  cf^  ■  cf^ 


Tlie  signal  received  by  a  station,  f ,  combines  with  a  station  transmitter  of  frequency,  f 
to  get  beat  frequency: 

dN 

—  =  f  -  f 
dt  » 


where  N  is  the  number  of  beats. 


1  .  1  .  1  . 
‘'s  -  =  ‘‘s  -  ‘ T  +  T  L,  +  -7  Lj  +—  L3 

®  '  c  ct.j.  ct  :j,  cf  ^ 


Define  Af  =  f  -  f_ 

s  T 


<m  ^  ^  ^  _L1]iL  I  I  clLa  ^  I  dLa 

dt  c  dt  cf.j.  dt  cf^  dt  cf^  dt 


1  1  1 

N  =  AfAt  +  -  +  —  AL,  +  —  ALj  -H—  ALj 

c  ct.^  ct  ^  cf  :j: 


Solving  for  the  range  difference,  Ap,  gives 


m 
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The  first  term  in  Equation  11,  represents  the  observed  uncorrected  value  for  the  range 

difference.  The  correction  terms  are  given  by  Api,  Ap2,  and  Apj. 


DUAL-FREQUENCY  DOPPLER  DATA 


Most  satellites  emit  two  signals  at  different  frequencies,  f j  ^  and  f^^.  to  eliminate  the 
first-order  correction. 


‘  T .  fi 


1 


fr  =  fr  -  -  -  “  L, 

c  cIt 

'  I 


L,  - 


1  ^.,-2  ^2  ,.f3  ^3 


Cf:^ 


I  .  p  >  • 

f  fr  =  f-r  -  fj  --  T~  L 

'2  ‘2  *  2  ,-r 


C  L  t  .j.  ^ 


1  •  I  • 

■  cri 

Tj  '2 


Assume  t'  <  f,.  and  define  X  =  f  /f  . 

‘1  ‘2  *1*2 


Subtracting  X  times  Equation  12  from  Equation  13  gives 


1  X\L, 


*r,  ~  ^*ri  ^  If 


1  X  \  L 


1  X  \  Lj 


‘I 


t  ^  I  ^  /  C 

T  T  ' 

‘  2  ‘  I 


Siibstitiitiiii:  for  X  eives 


t'  f. 


=  (f, 


-  f2 


't  /‘t 
1  2 


-  -  —1  p/c 

r. 


1  I 


4 


i  ,  'i  ,  'r,  J  c 


-3  ,-2 


'r.'rj 


Define  f  =  f 


>,  T, 


(12) 


(13) 


9 


y-  > . 
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Then  the  frequency  received  by  the  station  is 


To  determine  the  range  difference  for  dual-frequency  data,  the  same  analysis  as  for  the  single¬ 
frequency  case  is  followed.  The  result  (see  the  appendix)  is  given  below; 


Ap 


=  — —  (N  -  AfAt) 


(f_  +  f^  )  (f_  f^  ) 

'  1  2  I  2 


AL, 


NUMERICAL  ANALYSIS 


A  numerical  analysis  was  performed  to  determine  the  magnitude  of  the  second-  and 
third-order  corrections  to  the  phase  path  and  also  to  the  range  difference. 

A  FORTRAN  program  was  written  to  perform  the  calculations  and  interfaced  with  another 
program  that  contained  the  ARL  model.  The  combined  program  uses  the  Earth-fixed  coordinates 
of  the  satellite  and  the  given  station  position  to  determine  the  slant  range,  p.  at  given  time 
intervals  for  a  pass  over  the  station.  The  various  ionospheric  parameters  needed  by  the  ARL 
model,  along  with  the  magnetic  field  components,  are  evaluated  at  the  point  on  p  that  is 
350  km  above  the  Earth.  This  particular  height  was  chosen  after  numerical  integrations  showed 
it  to  give  the  best  average  fit  for  the  parameters  along  p.  A  simplified  fiowchart  of  the  com¬ 
bined  program  is  given  in  Figure  3. 

A  NAVSAT  orbit  provided  the  data  for  the  analysis.  Results  were  obtained  for  three 
stations: 

1.  Near  the  North  Pole 

2.  Mid-latitude 


ARE  COORDINATES  WITHIN 
AN  85  deg  ZENITH  ANGLE  OF 
STATION? 


YES 


COMPUTE  p 


t 


COMPUTER  IONOSPHERE 
COEFFICIENTS 


PERFORM  CLOSED  FORM 
INTEGRATIONS  FOR  CORRECTIONS 
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Station  Number 


Location 


118 

30280 

30121 


Thule,  Greenland 
Santiago,  Chile 
Quito,  Ecuador 


For  the  analysis,  a  relatively  high  sunspot  number  of  137  was  assumed.  For  the  single¬ 
frequency  case,  the  frequency  was  taken  to  be  150  MHz;  for  the  dual-frequency  case,  the 
assigned  frequencies  were  150  and  400  MHz. 


CORRECTIONS  TO  THE  PHASE  PATH 


For  single-frequency  data,  the  sum  of  the  second-  and  third-order  corrections  to  the 
phase  path  varied  from  -24  cm  to  -10  m,  depending  on  station  location  and  satellite  zenith 
angle.  For  two  frequencies  the  sum  of  the  corrections  is  less,  varying  from  -2.5  cm  to  -1.3  m. 

fhe  second-order  correction  was  found  to  be  much  larger  than  the  third-order  correction 
for  man>  points  along  the  pass.  This  was  especially  evident  for  station  118  (Thule.  Greenland) 
where,  for  small  satellite  zenith  angles,  the  second-order  correction  was  nearly  27  times  larger 
than  the  third-order  correction  for  the  single-frequency  case  and  almost  20  times  larger  for 
the  dual-frequency  case.  This  is  not  unexpected  since  the  magnetic  field  contribution  to  the 
precession  frequency  along  p  is  larger  here. 

For  the  single-frequency  case  at  station  30280  (mid-latitude),  the  second-order  correction 
was  approximately  5  times  larger  than  the  third-order  correction  for  many  points  along  the 
pass,  decreasing  to  1.5  times  as  large  at  other  points,  and  becoming  3.7  times  less  at  one 
point.  The  dual-frequency  case  showed  the  second-  and  third-order  corrections  to  be  about 
equal  in  magnitude  at  many  points  along  the  pass,  with  the  second-order  correction  becoming 
approximately  3.5  times  larger  than  the  third-order  correction  at  other  points,  and  5  times 
less  at  one  point  along  the  pass. 

The  results  for  station  30121  (equator)  showed  the  second-order  correction  to  be  from 
3.5  to  10.5  times  larger  than  tlie  third-order  correction  for  most  points  along  the  pass  and 
approximately  3  times  less  at  one  point  for  the  single-frequency  case.  For  two  frequencies,  the 
second-order  varied  from  2.5  to  7.6  times  larger  and  approximately  4  times  less  at  one  point. 

Tables  I  through  3  give  a  complete  listing  of  the  results  for  the  phase  path  corrections. 
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TABLE  1.  PHASE  PATH  CORRECTIONS. 
STATION  118 


Single  Frequency  Dual  Frequency 


Zenith  Angle 

L,  (ni) 

L2  (m) 

r 

i 

Lj  (m) 

Lj  (m) 

83.81 

-556.5431 

-1.5854 

-0.2581 

-0.1621 

-0.0363 

79.36 

-537.8851 

-1.7190 

-0.2290 

-0.1758 

-0.0322 

74.03 

-552.5074 

-2.0435 

-0.2312 

-0.2090 

-0.0325 

67.34 

-5 14.4646 

-2.2670 

-0.1976 

-0.2318 

-0.0278 

58.50 

-444.0260 

-2.3880 

-0.1548 

-0.2442 

-0.0218 

46.26 

-369.1418 

-2.4530 

-0.1208 

-0.2509 

-0.0170 

29.34 

-308.2193 

-2.4730 

-0.0999 

-0.2530 

-0.0140 

13.07 

-279.5940 

-2.4470 

-0.0908 

-0.2503 

-0.0128 

25.42 

-292.3383 

-2.3550 

-0.0915 

-0,2409 

-0.0129 

43.47 

-337.6948 

-2.2080 

-0.1022 

-0.2259 

-0.0144 

56.78 

-399.6440 

-2.0367 

-0.1242 

-0.2083 

-0.0175 

66.32 

-470.9486 

-1.8750 

-0.1609 

-0.1918 

-0.0226 

73.45 

-538.8791 

-1.7042 

-0,21  1  1 

-0,1743 

-0.0297 

79.07 

-594.6678 

-1.5250 

0.2710 

-0.1559 

-0.0381 

83.72 

-634.7061 

-1.3470 

-0.3304 

-0.1378 

-0.0465 
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TABLE  2.  PHASE  PATH  CORRECTIONS, 

STATION  30280 

Single  Frequency 

Dual  Frequency 

Zenith  Angle 

L,  (m)  Li  (m)  Lj  (m) 

Li  (m)  Lj  ( 

82.65 

78.37 

73.41 
67.49 
60.20 
50.94 
38.97 
23.99 

10.47 

19.62 

35.38 

48.42 

58.52 
66.41 
72.75 
78.01 

82.52 


-971.2966 

-1027.0405 

-1044.3015 

-1016.6601 

-953.4671 

-874.2020 

-793.8011 

-734.1684 

-725.51  16 
-788.6913 
-934.0459 
-1148.6793 
-1413.5384 
-1698.1758 
-1946.6436 
-2068.1360 
-1967.8629 


-4.3526 

-4.5764 

-4.6496 

-4.5264 

-4.2193 

-3.7739 

-3.1973 

-2.5000 

-1.6963 
-0.7774 
-0.2675 
- 1 .4604 
-2.7990 
-4.2300 
-5.5975 
-6.5857 
-6.7961 


-1.0043 
- 1 .0634 
-1.0451 
-0.9639 
-0.8585 
-0.7564 
-0.6733 
-0.626C 

-0.6399 

-0.7458 

-0.9817 

-1.3891 

-2.0193 

-2.8959 

-3.8845 

-4.5026 

-4.1003 


4452 

4680 

4755 

4629 

4315 

3860 

3270 

.2556 

1735 
0795 
0274 
1494 
2862 
4326 
,5725 
0.6735 
.(>95 1 


83.97 

-1040.172 

-6.08' 

79.31 

-91 1.7480 

-5.29 

73.71 

-779.4886 

-4.47. 

66.70 

-654.9483 

-3.68- 

57.54 

-544.7377 

-2.93. 

45.23 

-451.9222 

-2.19i 

29.65 

-386.7652 

- 1 .45- 

>7 

-1.1466 

0 

-0.8206 

19 

-0.5646 

•5 

-0.3920 

;i 

-0.2804 

>8 

-0.2089 

•3 

-0.1673 

13 

-0.1562 

18 

-0.1791 

>9 

-0.2386 

'5 

-0.3366 
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CORRECTIONS  TO  THE  RANGE  DIFFERENCE 


The  second-  and  third-order  range  difference  corrections  were  determined  by  taking  the 
difference  in  the  phase  path  corrections  at  two  successive  points  along  the  pass.  The  sum  of 
the  second-  and  third-order  range  difference  corrections  varied  from  as  little  as  0.05  cm  to 
- 1 .6  m  for  the  single-frequency  case  and  from  0.09  cm  to  approximately  30  cm  for  dual¬ 
frequency  data.  Tables  4  through  6  give  the  corrections  to  the  range  difference  for  the  three 
stations. 


TABLE  4.  RANGE  DIFFERENCE  RESULTS. 
STATION  118 


Zenith  Angles  for 

Which  Difference  Single  Frequency  Dual  Frequency 


was  Computed 

Api  (m) 

Ap2  (m) 

Apalm) 

Apj 

(m) 

Apj  (m) 

83.81, 

79.36 

18.6580 

-0.1337 

0.0291 

-0.0137 

0.4094 

X 

10-2 

79.36. 

74.03 

-14.6222 

-0.3244 

0.2180 

X 

10-2 

-0.0332 

-0.3066 

X 

10-2 

74.03. 

67.34 

38.0428 

-0.2234 

0.0336 

-0.0228 

0.4724 

X 

10-2 

67.34. 

58.50 

70.4386 

-0.1212 

0.0428 

-0.0124 

0.6020 

X 

10-2 

58.50. 

46.26 

74.8842 

-0.0651 

0.0339 

-0.6662 

X 

10-2 

0.4772 

X 

10-2 

46.26. 

29.34 

60.9225 

-0.0205 

0.0210 

-0.2101 

X 

10-2 

0.2956 

X 

10-2 

29.34. 

13.07 

28.6253 

0.0259 

0.9014 

X 

10-2 

0.2651 

X 

10-2 

0.1267 

X 

10-2 

13.07. 

25.43 

-  12.7442 

0.0925 

-0.6906 

X 

10-2 

0.9464 

X 

10-2 

-0.9713 

X 

10- 

25.43. 

43.47 

-45.3565 

0.1470 

-0.0106 

0.01503 

-0.1496 

X 

10-2 

43.47. 

56.78 

-61.9492 

0.1717 

-0.0221 

0.0176 

-0.3102 

X 

10-2 

56.78. 

66.32 

-71.3046 

0.1615 

-0.0367 

0.0165 

-0.5159 

X 

10-2 

66.32. 

73.45 

-67.9305 

0.1710 

-0.0501 

0.0175 

-0.7050 

X 

10-2 

73.45. 

79.07 

-55.7087 

0.1796 

-0.0599 

0.0184 

-0.8423 

X 

10-2 

79.07. 

83.72 

-40.0383 

0.1772 

-0.0594 

0.0181 

-0.8358 

X 

10-2 

Zenith  Angles 
for  Which  Difference 
was  Computed 

82.65,  78.37 

78.37,  73.41 

73.41.  67.49 
67.49,  60.20 
60.20,  50.94 
50.94,  38.97 
38.97,  23.99 
23.99.  10.47 
10.47,  19.62 
19.62,  35.38 

35.38.  48.42 

48.42.  58.52 
58.52.  66.41 
66.41.  72.75 
72.75.  78.01 
78.01.  82.52 
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TABLE  5.  RANGE  DIFFERENCE  RESULTS. 
STATION  30280 


Single  Frequency 


Dual  Frequency 


r -if.vjrf'.  ifi 


Api  (m) 

Ap2  (m) 

Apa  (m) 

Apz  (m) 

Apa  (m) 

55.7439 

-0.2238 

-0.0591 

-0.0229 

-0.8314  X 

10-2 

17.2610 

-0.0732 

0.0183 

-0.7489  X  10-2 

0.2577  X 

10-2 

27.6414 

0.1232 

0.0812 

0.0126 

0.0114 

63.1930 

0.3071 

0.1054 

0.0314 

0.0148 

79.2651 

0.4454 

0.1021 

0.0456 

0.0144 

■  *  ■ .  •  1  ■, 

80.4009 

0.5766 

0.0831 

0.0590 

0.0117 

59.6327 

0.6977 

0.0473 

0.0714 

0.6650  X 

10-2 

:vy.: 

8.6568 

0.8033 

-0.0139 

0.0822 

-0.1958  X 

10-2 

63.1797 

0.9189 

-0.1058 

0.0940 

-0.0149 

T"""m 

145.3545 

0.5098 

-0.2360 

0.0521 

-0.0332 

214.6334 

-1.193 

-0.4073 

-0.1220 

-0.0573 

264.8591 

-1.338 

-0.6302 

-0.1369 

-0.0886 

.  •/  %  %  •- 

284.6374 

-  1 .43 1 5 

-0.8766 

-0.1464 

-0.1233 

248.4678 

-1.3675 

-0.9886 

-0.1399 

-0.1390 

121.4924 

-0.9882 

-0.6181 

-0.101 1 

-0.0869 

100.2731 

-0.2103 

0.4023 

-0.0215 

0.0566 
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TABLE  6.  RANGE  DIFFERENCE 

STATION  30121 

RESULTS, 

Zenith  Angles 

for  Which  Difference 

Single  Frequency 

Dual  Frequency 

was  Computed 

Api  (m)  Ap2  (m)  Apj  (m) 

Ap2  (m)  Apj  (m) 

:■  83.97,  79.31 

1 28.4244 

0.7987 

:  79.31,  73.71 

132.2593 

0.8181 

73.71,  66.70 

124.5403 

0.7884 

1  66.70,  57.54 

1 10.2106 

0.7494 

57.54,  45.23 

92.8155 

0.7383 

:  45.23,  29.65 

65.1570 

0.7426 

1  29.65.  18.93 

18.1603 

0.7460 

1  18.93,  29.14 

-  35.2238 

0.6455 

1  29.14,  44.21 

-  81.5605 

-0.7971 

:  44.21.  56.13 

-111.5914 

-0.8175 

56.13.  65.03 

-125.1545 

-0.8124 

65.03.  71.86 

-137.0061 

-0.8117 

I  71.86.  77.33 

-159.8525 

-0.8651 

77.33.  81.90 

-207.4769 

- 1 .0343 

0.3166 


0.0458 

0.0837 

0.0360 

0.0806 

0.0243 

0.0766 

0.0157 

0.0755 

0.0101 

0.0759 

0.5844 

X 

10-2 

0.0763 

0.1564 

X 

10-2 

0.0660 

-0.3222 

X 

10-2 

0.0815 

-0.8364 

X 

10-2 

0.0836 

-0.0138 

0.0831 

-0.0199 

0.0830 

-0.0285 

0.0885 

-0.0445 

•0.1057 

-0.0770 

tz 


NSWC  TR  85-67 


CONCLUSIONS 


As  a  result  of  this  study,  the  authors  have  concluded  that  when  considering  higher  order 
ionospheric  effects,  both  second  and  third  order  must  be  considered;  neither  one  can  be 
neglected. 

In  the  near  future  the  authors  are  planning  to  merge  their  coding  for  the  ionosphere 
model  with  a  new  orbit  determination  program  (OMNIS)  that  is  currently  being  developed  by 
the  Space  and  Surface  Systems  Division  at  NSWC.  With  this  coding  included,  one  will  have 
the  option  of  computing  higher  order  effects  if  desired.  Also  in  the  rare  case  when  a  satellite 
loses  one  of  its  signals,  the  model  can  be  used  to  calculate  a  first-order  ionospheric  correction. 

A  technical  report  fully  describing  the  model  and  range  corrections  will  be  published  in 
the  near  future. 
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ALj  coefficient: 


ALj  coefficient: 


_ 1 _ 

fr,) 


Therefore 


Ap  =  -^  (N  -  AfAt)  + 


1 


AL,  + 
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